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Abstract 

Every normal toric ideal of codimension two is minimally generated by a 
Grobner basis with squarefree initial monomials. A polynomial time algorithm 
is presented for checking whether a toric ideal of fixed codimension is normal. 

1 Introduction 

Let ^ be a nonnegative integer d x n matrix of rank d, C{A) the cone in M'^ 
spanned by the columns of A, and the toric ideal associated to A as in [13] . 
The codimension of is n — d. If the columns of A form a Hilbert basis, 
i.e. C{A) n 7jA = f^A, then A and are called normal. Our main result is 

Theorem 1 Every normal toric ideal I_a of codimension 2 has a squarefree ini- 
tial ideal, and the corresponding reduced Grobner basis minimally generates I^. 

In Section 2 we introduce the necessary tools and prove Theorem [T] in the 
case of complete intersections. In Section 3 we complete the proof in the case 
of codimension two normal toric ideals which are not complete intersections. 
Section 4 deals with detecting normality and we prove the following result. 

Theorem 2 If codim(/_4) = n—d is fixed, then there is an algorithm to decide 
whether A is normal whose running time is polynomial in n and the bitsize of A. 

In other words, if n — d is fixed then one can decide in polynomial time 
whether n vectors in Z"^ form a Hilbert basis of the cone these vectors gener- 
ate. This answers a question raised by Alexander Barvinok at Snowbird (June 
2006) . We recently found out that Theorem [2] has been proved independently 
by F. Eisenbrand, A. Sebo, and G. Shmonin p]. The authors had announced 
the result at the 12th Combinatorial Optimization Workshop held in Aussois, 
France in January 2008. For recent advances on related complexity questions 
concerning lattice points in polyhedra of fixed dimension or codimension we 
refer to the article [4]. 

This work is motivated both by questions that are intrinsic to combinatorial 
commutative algebra and by applications to statistics and optimization. In the 
former domain, a longstanding conjecture states that every Cohen-Macaulay 
toric ideal has a monomial initial ideal that is also Cohen-Macaulay. This 
holds for toric ideals up to dimension 3 (see 0), and O'Shea and Thomas 



proved it for A-normal configurations [lOj. Theorem [T] offers supporting evi- 
dence because normal toric ideals and their square-free initial ideals are both 
Cohen-Macaulay. 

In algebraic statistics, our result ensures that the sequential importance sam- 
pling scheme of Chen, Dinwoodie and Sullivant [2] is applicable to exponential 
families with few states. In integer programming, it ensures that, for suitably 
chosen cost functions, every matrix of corank 2 specifies a Gomory family [7]. 
Finally, in algebraic geometry, where the definition of toric varieties requires 
them to be normal. Theorem [T] states that every toric variety of codimension 2 
admits a Grobner degeneration to a reduced union of coordinate subspaces. 

2 Complete Intersections 

The aim of this section is to prove Theorem [1] in the case when /_4 is a com- 
plete intersection. Before we approach this proof we will prove an alternative 
characterization of normality due to Sebo [12]. We include a proof of this result. 

Proposition 3 An integer matrix A = [oi, . . . , a„] is normal if and only if for 
each X G ker(^) there exists an integer vector y E ker(^) such that y < [x] . 

Proof. (=^) Suppose A is normal, and x G ker A. We have xioi + • • • -|-x„a„ = 
0. The vector z = \xi \ -01+ • • • + \xn] ■ an lies in the lattice ZA, and as 
z = ([xi] — xi) • ai + • • • -|- ([x„] — Xn) - an, it also lies in the cone C{A). 
Since A is normal, we conclude that z is in the semigroup N>1. We can write 
z = mitti -!-•■■+ mnCLn with mi, . . . , m„ nonnegative integers. The vector y 
with coordinates yi = \xi \ — mi lies in ker(^) and satisfies y < [x] . 

(<;^) Now suppose that for each x € ker(^) there is an integral y G ker(^) 
with y < [x] . Let z G C{A) H U^. This means that z = rioi -|- . . . -|- r„a„ with 
Ti G M>o and z = mioi + . . . + mnO-n with rrii G Z. Combining these we obtain: 

(mi - ri) • ai + • • • + {rrin - rn) • = 0. 

By hypothesis we may pick an integral y G ker(^) with y^ < \mi — r^] Then 

z = (mi - yi) • ai + • • • + (m„ - y„) • 

gives a nonnegative integral representation of z in terms of columns of A, since 
— yi > Q for all i. We conclude that z G N^, and hence A is normal. □ 

Our main tool in what follows is the Gale diagram of a vector configuration. 
Let ^ be a integer d x n matrix whose column vectors span M'^. We choose a 
matrix B whose rows form a lattice basis of ker(^) n Z". The set of column 
vectors of B is said to be a Gale diagram [15j of A. Normality of I_a is encoded 
in both A and in ker(^) = im{B), by Proposition [Sj and hence also in B. 

Hochster [6] proved that a normal toric ideal is Cohen-Macaulay. Thus 
a codimension two normal toric ideal has a minimal free resolution of length 
two. Peeva and Sturmfels |llj characterized Cohen-Macaulay codimension two 
lattice ideals. In this paper we consider saturated lattices whose lattice ideal is 
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Figure 1: The imbalanced Gale diagram of a complete intersection. 

a toric ideal. For the remainder of this paper we assume that is normal, or 
equivalently that ^ is a Hilbert basis of the cone C {A) . We assume that the 
cone C{A) is pointed and hence is homogeneous in some positive grading. 

The following result gives a supply of squarefree monomial terms of the 
binomial generators of This result has been proven in [13^, Proposition 4.1] 
and [9l Lemma 6.1], and we have also learned it from Winfried Bruns [Ij. 

Proposition 4 Suppose A is normal. Then each minimal binomial generator 
of the toric ideal I_a has at least one squarefree term. 

This implies that the conclusion of Theorem [1] holds when codim(/_4) = 1. 
In that case, I a is a principal ideal and the unique binomial generator of is 
a Grobner basis with its squarefree term being the leading monomial. 

In view of Proposition [4] our approach is to show the existence of a term 
order selecting the squarefree terms as initial terms. The Gale diagram gives 
information toward this goal. The following result is [11^ Proposition 4.1]. 

Proposition 5 // codim(/_4) = 2 then the following are equivalent: 

(i) The toric ideal Ij[ is not Cohen- Macaulay. 

(a) The toric ideal I_a has at least four minimal generators. 

(Hi) The matrix A has a Gale diagram B which intersects each of the four open 

quadrants in M^. Here the matrix B is identified with its set of columns. 

We now assume that is normal of codimension 2 and B is any Gale 
diagram. Then B = {{Bij,B2j) : j = 1, . . . ,n} intersects at most three open 
quadrants, and that any minimal generating set of I_4 has two or three elements. 
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In this section we examine the first case, where 1^^ = {x^ — x"^, — x^) is a 
complete intersection, and B is the 2 x n matrix whose rows are p — q and r — s. 
The Gale diagram B is called imhalanced if either Bij = or B2j > for all j. 

Lemma 6 \\\\ Theorem 3.9] A codimension 2 toric ideal Ia is a complete 
intersection if and only if there exists an imhalanced Gale diagram B. 

In the light of this lemma, we can represent a complete intersection by an 
imhalanced Gale diagram as depicted in Figure [TJ The arrows represent a sign 
class of columns of B and not just an individual vector. For instance the class 
labeled D in Figure [1] consists of all column vectors with Bij < and B2j > 0. 

Proof of Theorem\^ for complete intersections: By Proposition H] both genera- 
tors gi = x^ — x'^ and g2 = — have a squarefree term. The class of vectors 
F must exist in the Gale diagram since otherwise C{A) would not be pointed. 

If the term corresponding to F is squarefree then we can use any term 
order so that is the initial term of 52 and the squarefree term of gi is its 
initial term. Then gi and 52 forms the Grobner basis of /yi because their initial 
terms are relatively prime. This Grobner basis is or can be made reduced. 

Suppose that x^ is not squarefree. Then —B2j = / > 2 for some j € F, and 
B2j = 1 for j G i? U C U Without loss of generality we assume that x^ is the 
squarefree term of gi, so that By = 1 for j € Au B . We choose representatives 
from the D and E classes, labeling them —d and — e where d,e > 1: 
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Note that — //2 < —1. By Proposition [3] there exists an integral v G ker(^) 
with V < [n]. This vector is an integral combination v = a{p — q) + f3{r — s). 

If vi = then a = and v must be a positive multiple of (r — s) to ensure 
that vp = —Pf < r~//2l < —1- This implies the contradiction vb = P > 0- 

Next suppose that vi < —1. Then a < —1 (considering the A component) 
and P > 1 (considering the F component). The D representative requires that 
d + 1 < —ad + P < \^^~\ ■ But this implies that d = 0, a contradiction. We 
conclude that the D class is not present in the Gale diagram. By rotating 
the diagram by 90 degrees counterclockwise we can assume that we have an 
imhalanced Gale diagram where the B class is missing. For the two new minimal 
generators x^ — x'^ and x*" — x'' we are either in the first case analyzed above 
(i.e. X* and x^ are squarefree and relatively prime) or in the second case where 
x^ and x'^ are squarefree and relatively prime, as these do not contain any B 
variable. In both cases the two generators form a squarefree Grobner basis. □ 
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Figure 2: Gale diagram of a CM but not complete intersection configuration 

C 



D 


B2 






((0,1) 




E ^ 


,(0,0) 




A 




(1,0) 









A 



F2 G 

Figure 3: A syzygy triangle as in [TT] 

3 Normal but not complete intersection 

We now assume that the toric ideal is not a complete intersection, but it 
is normal and hence Cohen-Macaulay. This time we can assume that the Gale 
diagram is of the form as in Figure [2j As in Section 2, the vectors in the diagram 
represent a sign class of vectors. The minimal free resolution of I a has the form 



^ i?^ 



R R/Ia 



where R = . . . , A matrix representing the map R^ i?^ in the 

resolution can be determined using the recipe in [TT]. More precisely, using the 
two syzygy triangles in Figure 3 and Figure 4 we find that this matrix equals 



■ ABi 
B2CD 
F1F2 



D*E*F^ 
BIB^ 
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In the above matrix the letters represent classes of variables corresponding 
to the classes of vectors in the Gale diagram. Products of letters correspond to 
monomials in these classes of variables. The letters with an asterisk correspond 
to the same class as those without an asterisk, but they might have different 
exponent vectors (since they come from different syzygy triangles). The first 
column of the matrix corresponds to the first triangle and the second column 
to the second triangle. Moreover, by the Hilbert-Burch Theorem, the three 
2x2- minors of this matrix are precisely the minimal generators of 

ABiF^G* - B2CDD*E*Fl (1) 
ABiBlB^ - D*E*FiF^F2 (2) 
BIB2BICD - FiF2F;G* (3) 
The next result completes the proof of Theorem [TJ 

Lemma 7 There exists a term order such that the binomial generators , ([2]) j dS]) 
of the toric ideal form a Grobner basis with squarefree initial monomials. 

Proof. We have a few cases to consider. First, either the D class exists or it 
does not. If it exists then the monomial DD* is not squarefree and hence the 
squarefree term of ([T]) is the first term. Note that the first terms of ([2]) and 
(l3|) cannot be squarefree simultaneously: if they were, the Bi and B2 vectors 
cannot be present in the Gale diagram, and this would be an imbalanced Gale 
diagram. Similarly, the second terms of these binomials cannot be squarefree 
simultaneously. This gives two cases to consider. In the first case we have 

AG* - B2GDD*E*Fi, AB*2 - D*E*FiFi, B2B2GD - F^G* . (4) 
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Here Bi and F2 are absent because otherwise BiB^ and F2F2 are not squarefree. 
If we choose a lexicographic term order where A > G > {B2,C, D, E, Fi} 
then the underhned terms are the leading terms in The S-pair S{1,2) = 
D*E*FiF^G* - B2B^CDD*E*F{ is reduced to zero by the third binomial, and 
the S-pair S'(l,3) = AB2BICD - B2CDD*E*FiF* is reduced to zero by the 
second binomial. The S-pair 5(2, 3) reduces to zero since the leading terms AB2 
and FiG* are relatively prime, and hence ([3D is a squarefree Grobner basis. 
In the second case, the minimal generators and their squarefree terms are 



ABiF^G* - CDD*E*, ABiBl - D*E*F2 , BlCD - F2F2G*. 



(5) 



The product of the three underlined terms is equal to the product of the three 
non-underlined terms. Hence no term order selects the underlined terms as 
leading terms. However, the squarefreeness of these three monomials implies 
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From the diagonal edge in the two syzygy triangles we see that Bi = F2 = 1. 
This means that the non-underlined terms of the second and third binomials in 
([5]) are actually squarefree, and we are back in the previous case (jl]). 

Now suppose that the D vectors are not present in the Gale diagram depicted 
in Figure [2j Then the binomial generators ([T|) , ^ and ([3]) have the form 

ABiF^G* - B2GE*F^, ABiB^B^ - E*FiF^F2, BIB2BIC - FiF2F^G*. 

If in the first binomial the first term is squarefree we are back to ([!]) or (l5|) . If the 
second term is squarefree, then we rotate the Gale diagram 180 degrees. This 
leads to the same binomials but now with the first term of the first binomial 
squarefree. Once again we are back to (|4]) or ([5]). This concludes the proof. □ 



4 Checking Normality 

In this section we assume that the codimension m = n — d of is fixed. First 
we reformulate Proposition [3l Let z be an integral vector in ker(^). We define 

Pz = { 2; G ker(^) : [xj] > Zi for i = 1, . . . , n}. 

Since Pz = {a; € ker(^) : Xi > Zi — 1 i = 1, . . . ,n} and since we assume that 
the cone C{A) is pointed, Pz is a relatively open polytope in ker(^) ~ M™. 

Remark 8 If u and z are lattice vectors in ker(^) then Pu+z = z + Pu- 

Now let ;B be an n X m matrix whose columns form a lattice basis of ker(^), 
and let 6, be the rows of B. Then Pz is affinely isomorphic to Qv = {y ^ 1^™" '■ 
bi ■ y > bi ■ V — 1, i = 1, . . . ,n} where v is the unique lattice point in Z*" such 
that Bv = z. Remark [8] implies that for two lattice points v and w in Z"^ we 
have Qv+w = 10 + Qy. Note that Qo = {y £ K"^ : bi ■ y > —1, i = 1, . . . , n}. 
We now see that the following is equivalent to Proposition [3l 
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Theorem 9 The toric ideal Ij( is normal if and only if Qq + Z"* = M"* . 

Given any polytope Q = {y £ M"^ : Cy > d} of dimension m the smallest 
positive real number t such that tQ + Z™ = is called the covering radius 
of Q. If Q is a rational polytope it is known that the covering radius of Q is a 
rational number with a bit-size that is a polynomial in the bit-size of C and d. 

Corollary 10 The toric ideal is normal if and only if the covering radius 
of Qo, the closure of the polytope Qq, is less than 1. 

Proof of Theorem \^ Ravi Kannan [51 Section 5] has shown that, for fixed m, 
and given a rational m-dimensional polytope Q = {y E : Cy > d} where 
C G Z"^™ and d (z Z^, there exists an algorithm to find the covering radius 
of Q with runtime a polynomial in n and the bit-size of C and the vector d. 
Since one can compute a B whose bit-size is a polynomial in the bit-size of A 
in polynomial time, the above corollary implies the result. □ 
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